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Введение. Предлагается новая схема плоской статически определимой регулярной решетки. Стержни решетки 
соединены шарнирами. Цель исследования — вывести формулу зависимости от числа панелей первой частоты 
собственных колебаний узлов, наделенных массами, каждая из которых имеет две степени свободы в плоскости 
решетки. Жесткость всех стержней принята одинаковой, опоры (подвижный и неподвижный шарниры) неде- 
формируемыми. Другая цель исследования — найти в аналитической форме зависимость усилий в наиболее 
сжатых и растянутых стержнях от числа панелей. 

Материалы и методы. Используется приближенный метод Донкерлея определения нижней оценки собствен- 
ной частоты колебаний решетки. Жесткость решетки находится в аналитической форме по формуле Максвел- 
ла — Мора. Усилия в стержнях и реакции опор определяются из уравнений равновесия, составленных для всех 
узлов решетки. Обобщение результата на произвольное число панелей выполняется методом индукции с при- 
менением операторов символьной математики Мар|е по аналитическим решениям ряда задач для решеток с 
различным числом панелей. 

Результаты исследования. Нижняя аналитическая оценка первой частоты колебаний хорошо согласуется с 
численным решением для минимальной частоты спектра колебаний конструкции. Найдены формулы для уси- 
лий в четырех наиболее сжатых и растянутых стержнях и их линейные асимптотики. Все необходимые преоб- 
разования производятся в системе символьной математики Мар(е. 

Обсуждение и заключения. Полученная зависимость первой частоты колебаний решетки от числа панелей, 
массы и размеров конструкции имеет компактную форму и может быть использована как тестовая задача для 


численных решений и при оптимизации конструкции. 
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Введение. Решетчатые конструкции широко применяются в машиностроении как несущие или ограждающие 
элементы. Методика расчета собственных частот колебаний решетки традиционно базируется на численных 
расчетах, как правило, использующих метод конечных элементов с применением специализированных пакетов. 
Значительно реже применяются аналитические методы. Развитие математических пакетов символьной матема- 
тики (Маре, \УоШат Мафетанса, Оепуе, Махипа и др.) дает возможность искать такие решения для регуляр- 
ных систем, в которых возможен учет порядка регулярности (числа периодических структур конструкции, 
например, числа панелей) в решении. Цель исследования — найти аналитическую зависимость нижней грани- 
цы первой частоты решетчатой конструкции от числа панелей. Зависимость решения от числа панелей значи- 
тельно расширяет область применения формулы и дает возможность оптимизировать конструкцию по какому- 
либо параметру. Впервые проблемой существования и расчета статически определимых регулярных стержне- 
вых систем занялись Хатчинсон Р. и Флек Н. [1, 2]. В [3, 4] такие конструкции изучались в связи с задачами 
оптимизации. Существуют также аналитические решения в виде конечных формул для прогибов регулярных 
плоских [5—7] и пространственных ферм [8]. В справочнике [9] приводятся формулы для прогибов и смещений 
опор плоских ферм и решеток с произвольным числом панелей. В [10-13] получены нижние оценки первой 
собственной частоты колебаний плоских ферм, найденные методом индукции для произвольного порядка регу- 
лярной конструкции. 

Есть другое направление аналитических исследований конструкций [14—16]. В этих работах решение нахо- 
дится в системе Маре в виде тригонометрических рядов. 

В настоящей работе используется метод индукции, заключающийся в обобщении ряда отдельных расчетов 
решеток с последовательно увеличивающимся числом панелей на случай произвольного числа панелей. Во всех 
преобразованиях используются операторы системы символьной математики Маре. Объектом исследования 


Кирсанов М. Н. Аналитическая оценка частоты собственных колебаний плоской решетки 


является новая схема регулярной статически определимой плоской решетки прямоугольного очертания в виде 
несимметричной фермы на двух опорах. Ставится задача вывести аналитическую зависимость основной часто- 
ты колебаний конструкции от числа панелей. Выведенная формула может быть использована в задачах оптими- 
зации и для оценки численных решений подобных конструкций большого порядка, для которых в численных 
расчетах возможны погрешности, связанные с накоплением ошибок округления. 

Материалы и методы. В ферме длиной (2п-Та и высотой 4й содержится 5(п +1) узлов, включая и три 
опорных узла (один для крепления левой опорной стойки, два — для правых стержней, моделирующих непо- 
движный шарнир). Число стержней, включая три опорных стержня, М =10и+4. Ферма статически определи- 
мая. При расчете частот колебаний конструкции предполагается, что масса фермы сосредоточена в узлах. 


=== 


11 2и+1 


б и1 х 


Рис. 2. Нумерация узлов решетки, и=5 (рисунок автора) 


Вычисление усилий выполняется в системе Маре по программе [17]. Узлы и стержни фермы нумеруют- 
ся (рис. 2). Начало координат находится в левой опоре. Координаты задаются в циклах. 

Структура решетки устанавливается порядком соединения стержней в узлах. Для этого вводятся специаль- 
ные списки Ф, =[1,5ь] номеров й,Ь концов стержней а=1,..,№М. Стержни нижнего внешнего контура, 
например, имеют следующие номера узлов по концам: Ф, =[1,1+1,1=1,.., п. Таким же образом задаются и но- 
мера концов остальных стержней решетки. Система уравнений равновесия узлов в проекциях на оси координат 
составляется в матричном виде @5 =В, 5$ — вектор всех усилий в стержнях, включая и три реакции опор, 
В — вектор внешних узловых нагрузок. Матрица системы С состоит из направляющих косинусов усилий. 
При этом одно и то же усилие приложено к разным концам стержня и направлено в разные стороны: 

тер = -: 11 Сю а = - 11; 
Сы 2 бе =, И, 1=Ъ...,М, 


где 1, = Хх №, ‚ | = Уф, — Ух. , — проекции условных векторов стержней на оси координат; 


ы 1,2 


Машиностроение и машиноведение 


а. +. — длина стержня #=1...М№. 


Усилия в стержнях решетки можно получить из решения системы уравнений в символьном или численном 


виде. 317 
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Результаты исследования. Рассмотрим напряженное состояние решетки в случае нагружения по всем уз- 
лам вертикальными силами Р (рис. 1). На рис. 3 дана картина распределения усилий в стержнях конструкции. 
Толщины отрезков на рисунке условно пропорциональны модулям соответствующих усилий. Синим цветом 
выделены сжатые стержни, красным — растянутые. Значения усилий отнесены к величине Р нагрузки на узел с 
округлением до двух значащих цифр. Наиболее растянутый стержень ожидаемо оказался в середине нижнего 


пояса, наиболее сжатый — в нижнем стержне на правой боковой стороне конструкции. 


-13 


3,3 7,1 0, 8,3 5,5 
Рис. 3. Распределение усилий по стержням решетки, а=4 м, 6=3 м, п=5 (рисунок автора) 


Аналитические зависимости усилий в наиболее сжатых и растянутых стержнях от числа панелей получают- 
ся методом индукции из обобщения последовательностей отдельных решений. Например, для усилия У, В 


нижнем стержне левой стороны решетки последовательность значений для решеток порядка и = 1, 2, 3, ... имеет 
вид: У,/Р=-2, -12, —28, —54, —84, -—126, —170, —228, —286, —360,.... Рекуррентное уравнение для об- 


1п 
щего члена этой последовательности дает оператор гоЁ Нпагесиг системы Маре: 
У =У 1 2 „_ 2 из У —4 | И „5. 


‚п 1,и— 


‚п 


Решение уравнения с помощью оператора гзо]уе: 


У, =-Р(30и? —2(7-(-1")п- (-1" +1 /(8@2п-1). 


1,п 


Таким же образом находятся и другие выражения для критических усилий: 
О, =-Ра(20п? —30п? +5—4(с0$ ф+зтф)-2п)—с0$(26)) / (32(2п—1) п); 
И, = —0О,; 


у,, =-Р(би? —п-1)/(2п-1), 


гдеф=ли/2. Интересно заметить, что усилия У, и \, „ не зависят от размеров а ий. 


Аналитическая форма решений позволяет найти их асимптотики с помощью оператора Пти системы Маре: 
Шт У, „ =-ЗР.„ /8, ШПУ, „ =-ЗР. / 8; 
п>о 


1,п 5ит 5ит 
п 


ЫпО, /п=—Ш0, /п=-аР.„ / (161), 


5ит 
по п—>о5 


где Р„ =5"Р — суммарная нагрузка на решетку. Для усилий \/,„ и \, „в стержнях на боковых сторонах ре- 


вит 
шетки асимптоты горизонтальные, для усилий О, и (, на верхнем и нижнем поясах — наклонные. 

Частота собственных колебаний. Из всего спектра частот собственных колебаний конструкции для оценки 
ее динамического поведения наиболее важной является первая, низшая, частота. Ее значение входит в боль- 
шинство решений задач динамики конструкции. Эта величина требуется в том числе и для оценки сейсмиче- 
ских характеристик сооружения. Нижнюю границу первой частоты для регулярных конструкций в форме зави- 
симости от числа панелей можно получить аналитически. 
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При определении собственных частот колебаний конструкции принята упрощенная модель инерционных 
свойств фермы. Предполагается, что стержни решетки не имеют массы, а вся масса распределена равномерно 
по узлам. Пренебрегая движением опор, получаем общее число степеней свободы равное К =10и. Формула 
Донкерлея [13] для оценки нижней границы первой частоты имеет вид: 


К 
=», (1) 
р= 
где &, — парциальные частоты конструкции. Парциальные частоты колебаний масс определяются из 
уравнения: 
ту, +Бьу, =0, р=1....К. (3) 
Здесь у, = у,() — координата узла р; у, — ускорение; ), — жесткость, величина обратная податливости 


6, =1/Б, . Податливость можно вычислить по формуле Максвелла-Мора: 


к (р) 2 
6, =1/2,=>.(5%) 1, (ЕЕ), (4) 
а=1 
где ты — усилие в стержне с номером © от действия вертикальной единичной силы, приложенной к узлу р, 


где расположена масса. Коэффициент жесткости и парциальная частота зависят от места, где расположена мас- 
са. Для гармонических колебаний у, =И, т(!+ф) из (3) следует, = |, / т. Подстановка этого выраже- 


ния в (4) дает формулу для оценки первой частоты только по парциальным частотам колебаний масс: 


К К 
вр => в =т>. 5, =иКА,,+А,,). (5) 


р= р=1 


Отдельно выделены суммы для колебаний по вертикальному направлению Д,, и горизонтальному Д,,. 
Последовательный расчет частот колебаний решеток различного порядка показывает, что коэффициент Д,„, 
в (5) имеет вид: 

А, , = (141) / (ЕЁ); 

А, = (44а? +75с* +2501) / (9ЕЕЙ”); 

Д., = (24а? +25с* +361) /(ЕЕЙ); 

А, , = 2154а* +175с° +176^) / (1ЕЕ/”); 

А; = 2(1254а? +1125с° +7911°) /(27ЕЕР). 


В общем случае для парциальных частот по вертикали: 


А, =(Са’+С,с'+С,’) (ЕЕ). (6) 


Для парциальных частот по горизонтали: 


А, = (Са? +С;с* +1?) /(ЕЕа?). @ 


Используя оператор гэЁ_Япагесиг системы Маре, получаем однородное рекуррентное уравнение седьмого 


порядка для определения коэффициента при а? : 


С, 


1,п 


=ЗС С „2 5С! п С па С „5 ЗС! „-в С „1 ы 


Решение уравнения дает оператор гзо]уе: 


С = (и —2п? +61? —6(3(—1)" +10)п +9(-1)" —9) /18. (8) 
Таким же образом находятся и другие коэффициенты: 


С, =25(п-1/6; 


9 
С; = 36п? + 2(9(—1)" -14)п-9(-1)" -15)/ (6(2п-1)). © 


Аналогично: 
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С. = (164и? — (36 с0$(26) + 142)? + (24с0$(26) —245тф-+ 24с0$ф-+ 20)и-— 
—3с0$(20) +12 5тф-—12с0$ф+ 27) / (2—1) /12; 
С; = 10° —5п? +151 +8) / (2(2п-—1); (10) 
Се = (36? — (28+8с05(2ф))п? + (128+22с03(2) + 24со5ф-245тф)п-— 
—12с0зф-9с03(26)+109+12 т ф) /(2и-—1? /2, 


где ф=ли/2. 
В результате из (5—7) следует выражение для нижней оценки первой частоты: 
р =т((Са* + С,с°+ С?) + (С° + С° + СЫР) 1а? (ЕЕ). (1) 


с коэффициентами (8), (9) и (10). 

Оценка погрешности решения (11) возможна из сравнения с минимальной частотой всего спектра собствен- 
ных частот решетки, полученной численно. Спектр системы с многими степенями свободы находится из реше- 
ния задачи о собственных числах матрицы. Дифференциальные уравнения динамики масс конструкции с чис- 


лом степеней свободы К записываются в матричном виде: 
ти У+Р.У=0 (12) 
где Ох — матрица жесткости фермы; У — вектор смещений масс; |, — единичная матрица. Пусть Вх — мат- 
рица, обратная О„. Умножение (12) слева на Вх, дает уравнение: 
тв, У+ЬУ=0. (13) 
Связь У =-®”`У справедлива для гармонических колебаний с частотой ®. Отсюда из (13) следует задача о 
собственных значениях: В„ У =ЛУ, где Л, =1/ (т) — собственные значения матрицы Ву. Элементы матрицы 


податливости находятся по формуле Максвелла-Мора. Для вычисления собственных значений матрицы в си- 
стеме Мар[е предназначен оператор Е1зепуаез из пакета линейной алгебры ГлпеагА1юебга. Решение этой зада- 
чи в общем случае можно получить только в численной форме. 


Пример. Стальная решетка имеет длину панели а = 2 м, высоту й =1м, массы в узлах т = 300 кг, жесткость 


хх 5 хх хх 
стержней ЕЁ = 0,8.10`кН . На рис. 4 показаны зависимости первой частоты от числа панелей, полученные чис- 


ленно и аналитически. 
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Рис. 4. Сравнение первой частоты собственных колебаний решетки ©, и ее нижней аналитической оценки @р 


в зависимости от числа панелей (рисунок автора) 


Кирсанов М. Н. Аналитическая оценка частоты собственных колебаний плоской решетки 


Относительная погрешность 8 = (® —@©,)/@,, в зависимости от числа панелей (рис. 5), показывает, что с 


увеличением числа панелей погрешность, начиная с п =1, растет, а затем монотонно и достаточно быстро убы- 
вает. Это особенно важно при использовании полученной формулы при расчете решеток с большим числом 
панелей, для которых численный счет начинает набирать погрешность, связанную с накоплением ошибок 
округления, а затраты на компьютерные ресурсы быстро растут. 


= 


0,35 


0,30 


0,25 


0,15 


1 2 3 4 5 6 7 8 9 01113 м 15 106 17 Вт 


Рис. 5. Погрешность аналитического решения для собственной частоты в зависимости от числа панелей и высоты й 
(рисунок автора) 

Обсуждение и заключения. Предложена схема статически определимой плоской решетки. Поставлена за- 
дача получить аналитическое выражение нижней границы первой частоты собственных колебаний решетчатой 
фермы для произвольного числа панелей в конструкции. Решение получено методом индукции в системе 
Маре. Наличие экстремумов на построенных кривых дает возможность оптимизировать число панелей решет- 
ки, выбрав наибольшую точность оценки и подобрав требуемую частоту колебаний. В предложенном исследо- 
вании учтены горизонтальные колебания масс. Учет горизонтальных колебаний несколько усложняет итоговую 
расчетную формулу, делая ее более громоздкой. Помимо использованного метода Донкерлея, для оценки пер- 
вой частоты есть более точный энергетический метод Рэлея, дающий оценку первой частоты сверху. Однако и 
это решение в рассматриваемом случае имеет излишне громоздкий вид и здесь не приводится. 
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